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Abstract 

A classical dynamical model that treats break-up stochastically is presented for low energy 
reactions of weakly-bound nuclei. The three-dimensional model allows a consistent calculation 
of breakup, incomplete and complete fusion cross sections. The model is assessed by comparing 
the breakup observables with CDCC quantum mechanical predictions, which are found to be in 
reasonable agreement. Through the model, it is demonstrated that the breakup probability of the 
projectile as a function of its distance from the target is of primary importance for understanding 
complete and incomplete fusion at energies near the Coulomb barrier. 

PACS numbers: 25.70.Jj, 25.70.Mn 
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Recent developments of radioactive isotope accelerators provide an opportunity to inves- 
tigate on Earth the fusion reactions that form heavy elements in the cosmos. These involve 
reactions of nuclei far from stability, the most exotic of which are often very weakly bound. 
Breakup of weakly bound nuclei is thus an important process in their collisions with other 
nuclei. A major consequence of breakup is that not all the resulting breakup fragments 
might be captured by the target, termed incomplete fusion (icf); capture of the entire pro- 
jectile by the target is called complete fusion (cf). Such icf processes can dramatically 
change the nature of the reaction products, as has been investigated in detail for the stable 
weakly-bound nuclei 9 Be and 6 ' 7 Li There, at energies above the fusion barrier, CF yields 
were found to be only ~ 2/3 of those expected, the remaining 1/3 being in icf products. 
Events where the projectile breaks up and none of the fragments are captured provide an 
important diagnostic of the reaction dynamics. This we call no-capture breakup (ncbu), 
also referred to as elastic breakup. 

In a conventional picture of fusion, two colliding nuclei will fuse if they overcome the 
potential barrier due to their mutual Coulomb and nuclear interactions. The additional 
breakup degrees of freedom when one of the colliding nuclei is weakly-bound makes the 
process very much more complicated. An outstanding theoretical challenge is to model 
the CF and ICF processes in such collisions, since this separation is crucial to understand 
the effects of breakup on fusion fl g Quantum nrechanica, few body approaches, such as 
the continuum discretised coupled-channels (CDCC) method [4|, |5[ and the time-dependent 



wave packet method (6(, cannot separate incomplete and complete fusion contributions to 
their absorptive cross sections [3] , since both result in depletion of the total few-body wave 
'unction. The CDCC method can, however, make reliable predictions of the ncbu process 
3], as will be exploited here. What then are the alternatives to the above models? A novel 
optical decoherence model has been suggested [3] but has yet to be implemented. Another 
approach is to use the concept of classical trajectories which allow CF and ICF events to be 
separated, as in the two-dimensional model of Ref. Classical dynamical models are also 
being successfully applied to describe atomic many-body collisions, including ionization Q], 
a process in some ways analogous to breakup. 

In this Letter we present a new three-dimensional classical dynamical reaction model that 
allows quantitative calculations of CF and icf yields above the barrier. This is achieved by 
the introduction of a stochastically sampled breakup function. Crucially, this function can 
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be determined from experimental sub-barrier breakup data [10|], or from CDCC calculations. 
We apply the model to two-body projectiles, although it can be extended to more complex 
many-body breakup scenarios. As individual trajectories are followed, the model allows a 
clear identification of CF, icf and ncbu events, allowing predictions of their yields as well 
as breakup fragment angular, kinetic energy and relative energy distributions. This unique 
capability provides a powerful tool for interpreting future fusion experiments with radioactive 
nuclei, and in applications to 7-ray spectroscopy The model is assessed quantitatively 
by comparing its predictions for ncbu observables with those of fully quantum mechanical 
CDCC calculations. These are carried out for a simplified test case. The classical model CF 
and ICF predictions are then compared to experimental data. 

The classical model calculations are as follows. The weakly bound (two-body) projectile 
P, with incident energy Eq and orbital angular momentum Lo, is incident on the target T, 
initially at rest in the laboratory frame. Prior to breakup the projectile follows a classical 
orbit with a P-T distance of closest approach R min (E , L ). These orbits of the bound 
projectile are calculated by numerical solution of the classical equations of motion in the 
presence of the Coulomb and nuclear potentials between P and T. These generate a Coulomb 
barrier for head-on (L = 0) collisions of height V£ T at a separation Rg T . 

The effects of the Coulomb and nuclear two-body interactions that cause breakup are 
encoded in a local breakup probability V^^R), a function of the projectile-target separa- 
tion R, such that V B ' u (R)dR is the probability of breakup on the interval R to R + dR. 
Crucially, for a given projectile-target combination, both experimental data [10] and CDCC 
calculations indicate that the integral of this breakup probability along a given classical 
orbit is an exponential function of its distance of closest approach, R min (E , L ). Explicitly, 

Vljj{R)dR = A exp(-aR min ), (1) 

from which it follows uniquely that the local function has the same exponential form, 
T>bu( r ) K exp(-a,R). 

For each L Q (chosen to be an integer number of H) a sample of N incident projectiles is 
taken. The position of breakup on this orbit is determined by sampling a breakup radius 
Rbu on the interval [R m i n (E , L ), 00] with the weighting Vbu( r )- Although dependent on 
the constants A and a, the exponential weighting of Vbu(R) will clearly place most Rbu 
in the vicinity of R m in- If the chosen L is less than the critical partial wave for projectile 
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fusion, L cr , then the associated trajectory would normally lead to CF, i.e. R m in < Rb ■ F° r 
these L we set R m i n = R^b \ when sampling Rbu, and all breakup events are confined to the 
incoming branch of the projectile trajectory. For L > L cr breakup can take place on both 
the entrance and the exit branches of the classical orbit, which are sampled equally. Having 
chosen the position with the projectile center of mass at breakup, it is instantaneously broken 
up into fragments Fl and F2. 

Following breakup, fragments Fl and F2 now interact with T, and with each other, 
via real central two-body potentials having Coulomb barriers V B 3 at separations Rg, i,j = 
1, 2, T, i 7^ j. The instantaneous dynamical variables of the excited projectile at breakup, its 
total internal energy e 12 , angular momentum £ 12 and the separation of the fragments rf 12 are 
all Monte Carlo sampled. The initial separations d± 2 between the fragments are Gaussian 
distributed in their classically allowed region to mimic the radial probability distribution 
of the projectile ground-state (g.s.) wave function. For high £ 12 excitations, when there 
is no barrier between Fl and F2, d\ 2 is equated with their external turning point. In the 

— * 

test calculation below, the orientation of d\ 2 is chosen randomly over the 47r solid angle, 

— * — # 

whilst the orientation of £\ 2 is chosen randomly from all directions orthogonal to d 12 . Other 
distributions could easily be encoded in cases where the two breakup fragments are not 
identical. £\ 2 is sampled uniformly on the interval [0,£ maa; ]. For E\ 2 we tested two sampling 
functions for energies between the top of the barrier, V^ 2 , and a chosen maximum e rnax . Using 
uniform weighting or an exponentially decreasing weighting yielded very similar outcomes, 
although the convergence of the observables was faster for the latter, which was therefore 
used. Both £ max and e max were increased until convergence of the observables occured, as 
in the CDCC calculations. 

Having fixed the positions and dynamical variables of the excited projectile fragments 
at the moment of breakup, the instantaneous velocities of the particles Fl, F2 and T are 
dictated by conservation of energy, linear momentum and angular momentum in the overall 
center of mass frame of the projectile and target system. These breakup initial conditions 
were transformed to the laboratory frame where the equations of motion are solved. The 
computed trajectories of Fl, F2 and T determine the number of ICF, CF and ncbu events, 
fragment Fj being assumed to be captured if the classical trajectories take it within the 
fragment-target barrier radius R^ . 

From the N breakup events sampled for each projectile angular momentum L , the 
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numbers of events Ni in which i — (ncbu), 1 (icf), or 2 (cf) fragments are captured 
determine the relative yields Pi = Ni/N of these three processes after breakup, with P + 
Pi + P 2 = 1. The absolute probabilities Pi(E ,L ) of the three reaction processes are 
determined from the relative yields and the integrated breakup probability over the whole 
trajectory P B u(Rmin)- 

Po(Eo,L ) = PBu(Rmin) Po, (2) 
Pl(E ,L ) = P B u(Rmin) Pi, (3) 
P2(Eq, Lq) = [1 — PBu(Rmin)] H(L cr — L ) 

+ PBu(Rmin) Pi, (4) 

where H(x), the Heaviside step function, is unity for x > 0. Finally, the physical cross 
sections are calculated using 

<Ti(Eo) = nX 2 ^(2L + l)Pi(E , L ) (5) 

Lo 

where A 2 = h 2 /[2mpEo] and mp is the projectile mass. Other observables, e.g. the angle, 
kinetic energy and relative energy distributions of the fragments from ncbu events, are 
calculated by tracking their trajectories to a large distance from the target (200 fm in the 
calculations presented below). 

In order to validate our approach, ncbu observables, reliably calculated within the CDCC 
quantum mechanical approach, are compared with the classical model calculation. We do 
so in a simplified model case for the reaction of a pseudo- 8 Be projectile P (assuming a 
weakly-bound state of two «-particles) with a 208 Pb target. The CDCC and classical model 
calculations are matched as closely as possible, as discussed below. 

We bind P in a nodeless s-state with a separation energy of 0.5 MeV in the Woods-Saxon 
(WS) potential well (V, r, a) = (16.7 MeV, 1.20 fm, 0.62 fm). The same potential was used 
to construct the continuum states of P, which are non-resonant in all partial waves. The 
Coulomb potential between all particles is that of uniformly charged spheres with a radius 
parameter 1.2 fm. The CDCC breakup model space includes (i) even £ 12 up to Ah, (ii) a 
maximum a-a relative energy of 12 MeV, divided into 12 bins in each partial wave, and (iii) 
ct-T potential multipoles up to quadrupole terms. The a-T optical potential is the real WS 
interaction (62 MeV, 1.39 fm, 0.62 fm). A diagonal absorptive WS potential (10 MeV, 1.30 
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FIG. 1: The Prtj yersus R m in from the CDCC model, the empirical line from the 9 Be+ 208 Pb 
sub-barrier data [1(3 ]. and the optimal function used in Fig. [3l 

fm, 0.74 fm) is included in all channels, its range chosen so that absorption occurs inside 
the fusion radius of the a-T system, i?^ T = 10.6 fm. Calculations were carried out at a 
laboratory energy of E = 45 MeV. 

A critical input to the classical model is the local projectile breakup probability V^u(R) 
at a separation R from the target, deduced from P BU (R min ). To allow comparison with the 
CDCC, PBui.Rm.in) has been deduced from the CDCC breakup cross sections, <tj(ncbu), 
for each P-T partial wave. This is well represented by an exponential function (beyond the 
strongly absorbed partial waves), with A = 1.38 and a = 0.293 fm -1 [see Eq. ([I])]. Experi- 
mental parameters defining PBu{Rmin) have been deduced from breakup measurements for 
the reaction 9 Be+ 208 Pb for a wide range of sub-barrier energies [h]]. The parameters which 
best describe the experimental ncbu yields are A = 1.68 x 10 4 and a = 0.922 fin -1 . The 
CDCC (dashed line) and empirical PBui.Rm.in) (thick solid line) distributions are compared 
in Fig. [TJ It is significant that the two functions are very different, the slope parameter 
a depending sensitively on the reaction and the breakup mechanism. For instance, our 
model of 8 Be breakup, unlike the physical 9 Be, does not generate Coulomb or nuclear dipole 
breakup couplings. Nor does it include the effect of neutron transfer channels, that may 
be significant in the empirical distribution 1^]. This comparison shows that the V^ v iR) 
have high diagnostic value, and many reaction observables are also sensitive to the assumed 
Vbu(R), as shown below. 

We can now use the CDCC-derived V^iR) as the sampling distribution in the stochastic 
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classical model calculations to assess the consistency of the observables calculated with the 
quantum mechanical and classical dynamical models. In sampling breakup l m ax = 4ft and 
£max = 6 MeV. Guided by the calculated g.s. probability of the projectile, the Gaussian 
distribution of d 12 values is taken to have a centroid and variance of 1.8 fm and 1.0 fm, 
respectively. The same real potentials are used as in the CDCC when solving the post- 
breakup three-body problem. The potential used for the P( 8 Be)-T system prior to breakup, 
(120.9 MeV, 1.39 fm, 0.76 fm), is consistent with the folding model potential derived from 
the a—T potentials and the g.s. wave function of P. In this way the extended size and weak 
binding of P are taken into account prior to breakup. Values of L < 100ft were found to 
be sufficient for the convergence of the observables. 

The calculated total absorption and total ncbu cross sections from the CDCC are 
cx(abs)=327 mb and <t(ncbu)=97 mb, compared with <r(abs)=300 mb and <t(ncbu)=93 
mb from the classical model. Fig. [2] compares the CDCC (thin solid curves) and the classi- 
cal model (thick solid curves) for (a) the breakup cross sections as a function of P-T partial 
wave, (b) the singles a-particle angular distribution, (c) the a-a relative energy distribution, 
and (d) the singles a-particle energy distribution. The classical and CDCC results are seen 
to agree quite well, both qualitatively and quantitatively, providing support for the classical 
dynamical treatment. 

Fig. [2] also shows the strong sensitivity of the (classical model) observables to the breakup 
function assumed. The dashed curves result from using the empirical Vbu(R) from the 
9 Be+ 208 Pb experiment, as opposed to the CDCC breakup function (thick solid curve). The 
former results in breakup occuring much nearer to R m in, and thus to the target. The ncbu 
spin distribution when using the empirical breakup function (panel a) now peaks at smaller 
angular momenta with a corresponding shift in the peak of the a angular distribution to 
larger angles (panel b). The strong tidal forces experienced by the non- fused a-particles then 
lead to an enhanced a-a yield of high relative energy (panel c) and corresponding changes 
to the singles a energies (panel d). 

These first results already indicate a potentially very useful, and very sensitive mapping 
between (the unobservable) Vbu(^) anc ^ the NCBU observables shown in Fig. [2j The breakup 
function, and reaction dynamics, probed by combining the classical model and experimental 
data, could thus provide much needed insight into the localisation and nature of the breakup 
and fusion mechanism for different reacting systems. 
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FIG. 2: The fully quantum mechanical CDCC model ncbu observables are shown by the thin 
solid curves. Those from the classical model using the same breakup probability function are quite 
similar (thick solid curves). Using an empirical V^u(R) deduced from the 9 Be+ 208 Pb reaction 
gives very different results (dashed curves), which are scaled by 0.1 where noted, demonstrating 
the sensitivity of the observables to V^jj(R). 

Finally, we apply these ideas to experimental data. Using Vbu(^) deduced from the 
empirical breakup function (thick solid line in Fig. [1]), the model predicts CF and ICF yields 
as shown by the thick solid and dashed curves in Fig. [31 These are close to the experimental 
data (triangles and circles). Only a small change in V^u(R) is required to give a very good 
reproduction of the data (thin lines in Figs. Q]and[5]). This result is the first quantitative 
result supporting the suggested linkage [Ifj] between the breakup process observed below 
the barrier and suppression of complete fusion observed at energies above the barrier. 

In summary, a three-dimensional classical dynamical model for low energy reactions of 
weakly-bound nuclei is presented, which for the first time allows consistent realistic calcula- 
tions relating below-barrier ncbu with above-barrier CF and ICF. In the test case studied, 
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FIG. 3: The classical model CF (dashed curves) and ICF (solid curves) cross sections for 8 Be+ 208 Pb 
are compared with the 9 Be+ 208 Pb experimental data of Ref. [1]. The energy and the cross section 
of the calculated and experimental data are normalized by the barrier energies (Bq = 40.7 MeV 
for thick curves, -Bo = 40.5 MeV for thin curves, and Bq = 37.5 MeV for experimental data) and 
their respective radii {Rq = 10.7 fm for theoretical curves, and R$ = 10.8 fm for experimental 
data), obtained by fitting the above-barrier total fusion (cf+icf) cross sections with the formula 
o fus = irRl(l - Bo/E c . m ). 

the model, which treats the break-up process stochastically, is shown to reproduce the differ- 
ential ncbu observables of the CDCC quantum approach quite accurately. We have shown 
that the breakup probability function V^ V {R) has high diagnostic value, being sensitive to 
the reaction mechanism and reactants, and that ncbu reaction observables are strongly sen- 
sitive to and can interrogate this model input. The model also demonstrates the importance 
of the breakup probability function, accessible from experiments at sub-barrier energies, to 
predictions of complete and incomplete fusion yields of weakly-bound projectiles at energies 
above the Coulomb barrier. 
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